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Integrable strongly correlated electron systems (see, e.g. ^) have been an important subject of research activity 
(see e.g. In a recent paper we proposed an integrable correlated electron model which is an eight-state 

version of the supersymmetric U model The latter has been extensively investigated by many authors Q. By 
eight-state we mean that at a given lattice site j of the length L there are eight possible electronic states: 

io), cijo), cyo), cyo), 

4i42lO)- sli^alO), 4243IO): 4ic],243l0), (1) 

where cj ^ {cj,a) denotes a fermionic creation (annihilation) operator which creates (annihilates) an electron of species 

a = 1, 2, 3 at site j; these operators satisfy the anti-commutation relations given by {c] ^,Cj^p} = SijSap- The 
eight-state supersymmetric U model contain one free parameter and has Lie superalgebra (//(Sjl) as its symmetry. 

In this letter, we promote the Lie superalgebra gl{'i\l) symmetry of the model to the quantum superalgebra 
C/q[(7/(3|l)] symmetry, thus giving a new integrable model of strongly correlated electrons with two free parame- 
ters. We then solve the two-parameter model by the coordinate space Bethe ansatz method and derive the Bethe 
ansatz equations. 

The Hamiltonian for our new model on a periodic lattice reads 

L 

H{g,K) ^iJjj+i(5,K), 

C / \ I sinhK / t t IN 
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2 cosh K sinh^ n / t t t 

(^cjjCj 2Cj^3Cj + l,3Cj-|-l,2Cj-|-l,l + h.C. 



sinh K{g + 1) sinh K(g + 2) 
Kn -Ko sinhK ^-^ 

2 cosh K(g -t- 1) sinh^ AC , , 
- sinh.(g+l)sinh.(g + 2) ("^'^"^-^"^'^ + -,+Mn,+i,2n,+i,3), (2) 



*E-mail: yzz@maths.uq.edu.au 
^E-mail: hqz@maths.uq.edu.au 



1 



where g, k are two free parameters, nj = rij^i + nj.2 + '^j,3 with nj,a = c-l^Cj^a being the number operator for the 
electron of species a at site j, 9{P — a) is a step function of {f3 — a) and 

^ = _ In ^i"^ ( = i In sinh^ «(g + 1) ^ = - In ^"^^ ^-9 (3) 

smhK{g + l)' 2 sinh k(7 sinh + 2) ' sinhK(5 + 2) 

The model contains correlated single-particle and pair hoppings, uncorrelated triple-particle hopping, and two- and 
three-particle on site interactions. As is seen below, the chemical potential terms are essential for the model to have 
the quantum superalgebra symmetry. It should be noted that for the periodic lattice chemical potentials can be 
dropped from the Hamiltonian without affecting the integrability. However, they play an essential role for an open 
lattice. As we have demonstrated in Q], without these terms, the boundary system would not be solvable by the 
coordinate Bethe ansatz method for a large class of integrable boundary conditions |^,|| . 

Some remarks are in order. For the choice of k = 0, the above model becomes the eight-state supersymmetric U 
model proposed in Q , whereas in the limit of 77 = k (first discarding the chemical potential terms in the Hamiltonian 
and then taking the limit), it reduces to the electron model introduced in 

J = l a [ J 

whose integrability is established in ]Tot . 

The model defined by (|^) is ?7g[(7^(3|l)] supersymmetric and is exactly solvable on the one-dimensional periodic 
lattice. This is because the local Hamiltonian Hj j^i{g, n) is actually derived through the quantum inverse scattering 
method using a [/q[(7^(3|l)] invariant R-matrix. To show this, we denote the generators of J7g[(7^(3|l)] by E^^, Mi ^ = 
1,2,3,4 with grading [1] — [2] = [3] = 0, [4] = 1. In a typical 8-dimensional representation V{K) of C/q[(7^(3|l)], 
the highest weight A — (0,0,015) itself of the representation depends on the free parameter g, thus giving rise to 
a one-parameter family of inequivalent irreps. Let {|a^)}^=i denote an orthonormal basis with |1), |5), |6), |7) even 
(bosonic) and |2), |3), |4), |8) odd (fermionic). Then the simple generators {-E/^}^=i and E'^'^'^Y^^^ are 8x8 

supermatrices of the form 

El = |3)(4| + |5)(6|, El = |4)(3| + |6)(5|, E\ = -|4)(4| - |6)(6| - |7)(7| - |8)(8|, 
El = |2)j3|+ |6)(7|, El = |3)(2| + |7)(6|, i?| = -|3)(3| - |5)(5| - |7)(7| - |8)(8|, 

[a], |1)(2| + ^[.9+l]J|3)(5| + |4)(6|) + ^[.9 + 2],|7)(8|, 
[g], |2)(1| + y/[.9+1]</(|5)(3| + |6)(4|) + ^[ff + 2], |8)(7|, 



El^ 
E:' 



^3 




£;| = -|2)(2|-|5)(5|-|6)(6|-|8)(8|, 

El^g\l){l\ + (g + 1) (|2)(2| + |3)(3||4)(4|) + (g + 2) (|5)(5| + |6)(6| + |7)(7|) + (g + 3) |8)(8|, (5) 

where [x]q ~ {q^ — <1^^)/{q — Q^^)- 

Uq{gl (3\1)) is a graded Hopf algebra with coproduct given by 

A{Ep^I^Ei: + Ei:^I, = 1,2,3,4, 

Associated with the 8-dimensional representation, there is a C/q[g^(3|l)]-invariant R-matrix which satisfies the graded 
Yang-Baxter equation. The R-matrix is given by 

Riu) - A + {2g)P2 + {2g){2g + 2)P3 + {2g){2g + 2){2g + 4)P4, (7) 

where (a) = (1 - g"+'')/(g" - q") and Pk : Pk[V{A) ® ^(A)] = V{Ak), k = 1,2,3,4, are four projection operators 
which now going to construct. Pi and P4 are 

Pi = j2\H){Hl P4 = El*t)(*fcl- (8) 

k=l fc=l 
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Throughout this paper, 

(*^| = , (la;) ® = {-l)i\-my)] (^y\ ^ (9) 

with [|a;)] = for even (bosonic) |a;) and [|a;)] = 1 for odd (fermionic) |a;). 1^^), I'^'l), k = 1, 2, ■ • • , 8 are given by 
I*1) = |1>®|1), 



(g^li)® i = 2,3,4, 



[^[g + W\^®\l)+q-'\l)®m 



|2) ®\i-2)-q-^i-2)^\2))l i = 5, 6, 

^J[9 + A,m^) ® |i) + i^'W ® |7)) + ^l\9\<,{q^\i) ® K) - <z-*|4) ® |3))], 

[\/[5]^(9-^+'|2) ® |7) + gi-i|7) ® |2) + gi+i|5> ® 4) 



+g-*-i|4) ® |5) - g-*|3) ® |6) - g*|6) ® |3)) + ^ [g + 2],{q^ \^) ® |1) + ® |8))], 



1 

yj{q3 +q-9)[2g+l]q 
1 



^(59+1 + g-s-l)(59+2 + g-9-2)[25 + 3]^ 

+q-i |5) |4) - qi |4) ® |5) + qi+^\2.) ® |6) - q-i-^\&} ® |3)) 

+V^(?^+'|1> ® |8) - «-^-'|8> ® |1))], 
1 



^[5 + 2],(5-*-2|7) ® |2) - gi+2|2) ^ |7) 



V(,.+2+,-.-^)[2, + 3]/-t^ + '^'^^' 
+-\/[5 + 2L(9^ |5) ® |6) - 16) ® |5))], 



-S-2| 



1 



-9-2| 



® |2) + g9+^|2)(g)|8)) 



(^|i)+gS+2|i)o|8)) 



1*^) = 



^{q9+^+q-a-^)[2g + Z]^' 
+05 + 2],(gi |z + 2) ® |7) - g-i |7) ® |i + 2))], i = 3, 4, 

^==l===(-g-«-^|8) |i) + g*+i|i> |8)), i = 5,6,7, 



(10) 



The construction of the projector P2 is a bit involved. We first construct the basis vectors for V{A2). With the 
help of the coproduct formulae and after some algebraic calculations, we get 

1 



^/¥Tq^ 
1 

y/q^Tq^ 

*3) = |2)®|2), 



(-g-*|2)0|l) + g*|l)®|2)), 



(-g-^|3)®|l)+g*|l)0|3)). 



(-g-5|4)0|l) + g5|i)^|4)), 



*^) = -^=i==(g-i |2) ® |3) + gi|3) ® |2)), 



^-e) = |3) ® |3), 



V9 + q-^ 
1 



(g-5|2)®|4)+g^|4)®|2)). 



(9-'|3)®|4)+g5|4)®|3)), 
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I*: 



^10 



*?1 



*?2 



13 



16 



*?8 



19 



22 



^23 



|4) ® |4), 



1 



1 

1 

1 
1 
1 



(-q-^-l|2)® 
(-q-*-l|3) ® 











1 






g-9)[2.gH 


-1]. 




1 






q-9)[2g^ 


-Ik 




1 





(-(7"2--L|2) ® 
[-q-i-^\i)® 



|5) + g^+i|5)® 
|5) + g*+i|5)® 
|6) + ® 
|6) + ® 
|7) + g*+i|7}® 
|7) + g*+i|7) ® 



|2)) 
|3)) 
|2)) 
|4)) 
|3)) 
|4)) 



|2) ® |3) + g-^-i |3) ® |2)) + ^/[.g + 1],(|1) ® |5) - |5) ® |1))], 
fc3+^|2) ® |4) + q-^-h\A) ® |2)) + + ® |6) - |6) ® |1))], 



V(g9 + g-9)[25 + l]g 
1 

1 

1 



W9^+i|3) ® |4) + g-^-i|4) ® |3)) + ^\a + ^,{\l) ® |7) - |7) ® |1))], 



= (g*+i|6) ® |3) - g"*"^|3) ® |6) |5) ® |4) - |4) |5))], 



(g5+'|7) ® |2) - g"*"'|2) ® |7) +g^|6) ® |3) - g-^|3} ® |6})], 



V (.-^-+.----)i2. + 3jJ V^(^-- 16) « 15) - gi 15) ® |6)) 
+ 0.g + 2],(g^+i|8) ® |2) + g-«-i|2) ® |8))], 
V(g^^^+g-'-^)[2. + 3]J V^(^-- 17) ^ |5) - g^ |5) « |7)) 
+ 05 + 2],(gS+i|8) ® |3) + g-«-i|3) ® |8))], 
V(g^^^+g-'-^)[2. + 3]J V^(^- 17) « |6) - g^ |6) « |7)) 
+ 05 + 2],(g^+^|8) ® |4) + g-3-i|4) ® |8))], 

^(gf+i|2)®|7)-g-*~i|7)0|2) 



^(g9 + g-9)(gf+i + g-f-i)[2.g + 1], 
+g^+i|5) ® |4) - g-^-i|4) ® |5) - g* |3) ® |6) - g'^ |6) ® |3)) 
+ 05 + 2],(-g*|8) ® |1) +g-*|l) ® |8))]. 



(11) 



However, this basis for ¥{^2) is not orthogonal, that is (^'fl^^) 7^ (5(m for some I, m. To make it orthogonal, we 
define a metric matrix gim by 



9lr. 



(vl/flvl/^), ;,m=l,2,-..,24, 



(12) 



where are defined by (|^). We then form a dual basis by means of the metric matrix, with basis vectors given by 



24 



m \ 1 



(13) 



where g'™ = {g ^)im is the inverse of the metric matrix gim- Then by construction, = Sim for all I, m, 

which implies that are orthogonal to l^*™) for all I, m. Thus the projector P2 is given by 



24 



A = E (14) 



k=l 



Finally the projector P3 is obtained through the relation Pi + _P2 + P3 + P4 = 1. 

On the L-fold tensor product space V ®V ® ■ ■ ■ ®V we denote R{u)jj+i = I'^U-i ^ ^(y) ^ j»(L-j-i)^ g^^^ define 
the local Hamiltonian by 



d 



u 



(15) 

tt=0 



We make the identifications: 

|1) = |0), |2) = c]j0), |3) = 4|0), |4)=43|0), 

|5)=c]i42l0), |6)-4i43|0), |7}=424,3l0), |8)=4i4243l0)- (16) 
Then by ( |l0| , ^1]) , (|9[ [l3|) and ([l^), and after tedious but straightforward manipulation, one gets, up to a 



constant. 



F„-+i(g,«) = = e-«). (17) 

K 



This identity also shows that H{g,K) commutes with the generators (^ of C/q[(7Z(3|l)], since the R-matrix R{u) is a 
C/q[(7^(3|l)] invariant. 

We now solve the system by means of the coordinate space Bethe ansatz technique. We assume the following 
wavefunction 

/ ^ \ 

(^Pqi > • • • , fcpQ„ )e-xp\i'Y^ kp^ Xj , (18) 

where Q is the permutation of the 'N particles such that 1 < xq^ < ■ ■ • < ^Qn ^ ^- Denote Xq = {xq-^ < ■ ■ • < ^Qn}- 
The coefficients AaQ_^ ,---,aQ^ (^Pqi i ' ' ' > ^Pqjv ) from regions other than Xq are connected with each other by elements 
of two-particle S-matrix: 

S'ij (/ci, fcj)°° = 1, 0=1,2,3, 

S^J{Kk,)t^-^^^^^^^^,. a b, a,b^ 1,2,3, 
sm(Ai Xj — IK) 

S.,ih,kXa - e-'g-t-^X^--^^) . . a,6 = 1,2,3, (19) 

sm(Ai — Aj — IK) 

where A^ are suitable particle rapidities related to the quasi-momenta kj of the electrons by 

k{\) = 2 arctan(cothctan A), (20) 

where the parameter c is defined by 

1 f, sinh i (77 + k) , 1 , , 

-i|-i .i..hi(„-») '-''l- 

The periodicity condition for the system on the finite interval (0, L) results in the Bethe equations for the set of N 
momenta kj : exp{ikjL) = Tj, j = 1, ■ ■ ■ , N , where 

^] = Sj.j+i{kj,kj+i) ■ ■ ■ Sj^N{kj,kN)Sj^i{kj,ki) ■ ■ ■ Sjj_i{kj,kj^i), j ^1,---,N. (22) 

The meaning of Tj is the scattering matrix of the j-th particle on the other (N — 1) particles. So now the problem is 
to diagonalize Tj to arrive at a system of scalar equations. It can be shown that Tj = t(A = kj), where 



5 



r(A) - tro [5o,i(A - h) ■ ■ ■ So,n{X - k^)] (23) 

is the transfer matrix of the inhomogeneous J7g[(7Z(3)]-spin magnet of N sites. The commutativity of the transfer matrix 
for different values of the spectral parameter A implies that Tj, j — 1, ■ ■ ■ , N can be diagonalized simultaneously. The 
Bcthe ansatz equations are written in terms of the rapidities A^j"^ , A^j"^ and A^. 



k,L ^ -Q sin(Aj- - Aj^' - iK/2) 



^J--^ sin(Aj - A^^^ + iK/2) ' 

yr sin(Ai,^^ - Aj + iK/2) _ sin(Ai^^ - A^^^ + ik) sin(Ai^^ - A^'^ - iK/2) ^ 

f}^ sin(Ai^) - A, - tK/2) ' y sin(A^^) - A^^^ _ ^^,) ^ sin(Ai^) - A^.'^ + zk/2) ' " ' ' " ' 
,(l)_L,■./o^ ^ o^.^a(2) _ a(2) 



^ sin(Ay^ - Ay^ + ^K/2) ^ sin(Ay^ - A^^> + ^K) 
y sin(A(^) - A« - Z./2) ysin(A(^)-A(^)-z.)' ^ " ^ ' 

The energy of the system in the state corresponding to the sets of solutions {A^^^ A^^-*} and {Ac} is (up to an additive 
constant, which we drop) E — — 2X]jLi cosfcj. 

To summarize, we have presented a new two-parameter integrable model which is an eight-state supersymmetric 
electron model with correlated single-particle and pair hoppings as well as uncorrelated triple-particle hopping. We 
have solved the model by the coordinate Bethe ansatz method and derived the Bethe ansatz equations. There are 
many problems remained to be done for this new model. One of them is to incorporate integrable boundary conditions 
into the model. We hope to report results on this aspect in future publications. 
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